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An  analysis  1 s  presented  of  the  aechanlsa*  it  host  transfer  and 
vaporisation  in  an  ldaalisad  two-dlaenslcnal  fual  apra/  consisting  of  two 
parallel  droplet  at  rasas  Injected  in  a  hot  gas*  With  the  asauaption  of  a 
constant  relative  velocity  between  the  gas  and  liquid  phases  (Inviscld 
approximation),  a  eye tea  of  partial  differential  equations  with  lines? 
operators  and  non-linear  source  tenia  Is  obtained*  With  the  use  of  Green's 
functions*  this  systaa  is  transforaed  Into  a  group  of  non-linear  ordinary 
differential  equations  and  integral  aquations  which  are  readily  solved  with 
auaerical  techniques.  The  diaenslonleas  groupe  influencing  the  problea  are 
Identified  and  a  discussion  is  presented  of  the  effects  of  four  relevant 
parameters  Including:  a  vaporisation-diffusion  tins  ratio*  a  liquid-gas  aass 
flow  ratio*  a  Spalding  transfer  nuaber*  and  a  Reynolds  nuaber.  ,  The  aethod 
yields  qualitatively  correct  eatlaatea  of  the  behavior  of  the  gaa  and  liquid 
phases  which  aey  serve  as  guidelines  in  the  developaent  of  costlier  nuaerlcal 
V  .detunes.  The  aethod  can  be  further  exploited  to  Include  chemical  reactions 
end  Co  ptedict  ignition  of  parallel  droplet  atreaas. 


Acftorislsn 
NI1S  '.lUt! 

DTK'  Tt* 

JuntKiceti  :n  _ 


By - 

Distribution/ 


Availability  Codes 
j Avail  nnd/o,* 
:iot  |  Special 


Introduction 


Tha  high  level  of  complexity  and  analytically  untractabla  non-linearitias 
present  in  tba  atudy  of  multidimensional,  ataady  and  unataady  spray  aodala 
bava  f oread  tha  combustion  scientist,  particularly  in  racant  years,  to  raly 
heavily  on  coatly,  larga-acala,  eoaputational  tachnlquaa  which  ara  abla  to 
handla  tha  rathav  challenging  aaount  of  mathematical  difficulties.  Examples 
of  thaaa  computational  afforta  have  appaarad  racantly  in  tha  literature*"*. 

To  aid  in  tha  understanding  of  tha  fundamentals  of  thaaa  phenomena,  it  ia 
appropriate  to  develop,  slmultanaoualy  with  the  numerical  schamae,  a  lass 
dostly  analytical  approach  based  on  simplified  models  that  retain  tha 
essential  physics  of  tha  problem.  A  basis. for  this  approach  has  bean  shown  in 
tha  linear  modal  of  Slrignano^. 

Xn  this  work  wa  establish  tha  fundamentals  of  a  multidimensional  spray 
modal  aimed  eventually  at  tha  study  of  tha  ignition  of  a  stream  of  droplets  by 
another  stream  which  has  already  ignited. 

He  consider  a  two-dimensional  uniform  velocity  gas  flow  in  which  tha 
vaporising  droplets  ar'  raprasanted  by  two  continuous  sheets  traveling  at  a 
different  velocity  within  tha  gas..  Uniform  heating  and  vaporisation  of  the 
droplets  is  considered  as  well  as  gas  diffusion  in  tha  transverse  direction. 
Axial  diffusion,  momentum  exchange,  density  variations  and  tha  affect  of 
convection  in  the  transverse  direction  are  neglected.  An  improvement  of  the 

model  to  account  for  spatial  variations  in  the  droplet  temperatures  cap  be 
achieved  with  the  use  of  eigen-function  expansions  end  integral  equations*" 

7.  The  chemical  reaction  terms  necessary  for  the  study  of  the  parallel-stream 
ignition  mechanisms  can  be  Incorporated  In  the  present  model  at  some  future 

point.  Realize  that  with  tha  use  of  Shvab-Zeldovich  variables  and  a  one-step 
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chemical  reaction,  all  but  one  of  tha  governing  equation*  la  the  reacting  caee 
..  can  be  caet  Into  a  non-raactive  form. 

Analysis 

la  thi*  work  we  consider  a  two-dimen* Iona  1 ,  cone  tent  property  gee  which 
flows  la  steady  state  between  two  planes  of  constant  but  different 
temperature.  Two  droplet  sheets  parallel  to  the  boundary  planes  are  flowing 
within  the  gas  with  the  saws  constant  velocity  0^  which  is  different  frow  the 
gas  velocity,  D(*  Tbs  fluid-wechanic  characteristics  of  the  problem  are 
simplified  by  assuming  that  the  velocities  and  are  constant  and  that  no 
anchange  of  momentum  takes  place  between  the  gas  the  the  liquid.  This  Is 
essentially  an  lnvieeld-flov  assumption  that  permits  the  solution  of  the 

•  .  s 

problem  In  terms  of  Green's  functions  and  Integral  equations. 

Density  variations  in  the  gas  phase  due  to  the  vaporisation  process  nay 
account  at  the  most  by  about  ten  percent  of  the  total  density  change.  Cas 
density  is  mainly  affected  by  temperature  changes  due  for  example  to  chemical 
reactions.  Accordingly,  neglecting  density  variations,  the  equations  of 
conservation  of  energy  and  conservation  of  the  fuel  species  ere: 
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where  a  and  D  are  the  thermal  and  mass  dlffusivities  respectively,  n  is  the 

number  of  droplets  per  unit  area  on  the  sheets,  p  is  the  gas  density 

and  4(x-Xj)  is  the  Dirac  delta  function  which  Indicates  that  the  source  teras 
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re  droplet  mass  an d  it  the  heat  of  vaporisation  of 

3  Is  art  defined  in  the  nomenclature • 
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Diaenalonless  variables  are  aalactad  aa: 
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Introducing  thaaa  dlaensionlasa  variablea  in  eqa.  (i)  and  (3)and  naglacting 
Cha  convactlva  tana  in  tha  tranavaraa  direction,  they  bacons: 


If  -F.  *Ji-  1:  v11*^  ,r-V>  v<rv 

at_  a2t.  ii 

ar-?i^r- 


(7a) 


(7b) 


ifl-  -  S 

d£  v 


(7c) 


<?d> 


where: 


S  •  R*  C_  in  (I  ♦  B 
v  Re  aff 


C  -  1  ♦  0.3  <K«)1/2  -  1  *  0.3  P4»')1/2 


Vr-T'«)  ■fps-?p 

*•«  *  s 


4 


The  boundary  condition!  become: 


At  C 

*  0: 

T*  -  T'^x), 

tf  *  V  r\  "  °*  v  ’  1 

(8a) 
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ar 
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At  X 

-  1: 

V  -  1, 

>yf 

»  o 

*X 

(8c) 

Equations  (7)  contain  tha  following  dimensionless  parameters  which 
characterize  tha  problea: 


(9) 

Other  paraaeters  are:  tha  relative  positions  of  tha  droplet  sheets,  the  ratio 
of  the  specific  heats  as  well  as  other  diaensionless  quantities  appearing 
through  the  equlllbriua  aquations: 
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(10b) 


We  concentrate  our  attention  on  the  effect  of  the  paraaeters  in  (9)  as  these 
are  very  influential  in  the  behavior  of  the  gas  and  liquid  phases  and  in 
particular,  in  the  development  of  the  fuel  concentration  field  in  the  gas 
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phase.  Th«  first  parameter.  Pj ,  can  b«  viewed  as  a  ratio  bstwaan  a 
characteristic  vaporization  tiae  and  a  characteristic  gas-diffusion  tiae.  If 
Pj  were  small ,  vaporization  of  the  droplets  would  take  place  in  a  tiae  whl:h 
is  shorter  than  ;  he  tiae  required  by  the  gas  diffusion  aechanisa  to  dlstrl  »ute 
the  fuel  over  the  cross  sectional  area.  In  the  Halt  Pj  ♦  o  there  will  be  two 
infinitely  narrow  bands  of  fuel  vapor  at  the  droplet  locations  and  essentlilly 
no  fuel  anywhere  else  in  the  gas.  As  Pj  increases,  diffusion  tends  to  smoith 
out  the  fuel  concentration  profiles.  The  second  paraaeter(  P2,  represents  the 
ratio  of  the  rates  of  incoming  liquid  aass  per  sheet  and  incoming  gas  mass  at 
the  channel  inlet.  P2  is  generally  lass  than  unity  and  as  it  increases,  tie 
result  is  an  Increase  in  the  local  and  total  aaount  of  fuel  in  the  gas 
phase.  The  third  paraaeter,  P-j ,  corresponds  to  the  Spalding  transfer  number 
based  on  the  reference  teaperatures  Tc  and  ,  T  ^  and  on  the  fuel  laf  nt  heat  of 
vaporization  L^.  The  last  parameter,  P4,  is  the  Reynolds  number  based  on 
initial  droplet  diameter  end  relative  velocity.  Solution  of  Eqs.  (7)  eubj set 
to  the  conditions  given  by  (8)  is  obtained  with  the  aid  of  the  Green's 
functions8  GT(£  ,x  |$ ' ,x ' )  and  CF<C»xlC'  *x'>  "h16*4  satisfy 

SG_  a2G 

jr  -  — r  -  «'<«  -  OS'Cx-x’) 


G^  •  0  at  x  •  0  and  x  •  1 


CT  •  0  for  C  <  ?' 


(11s) 
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Pi  — -Os’Cx-x’) 


3G?/3x  *  0  at  x  “0  and  x  "  1 


Cr  -  0  for  c  <  C'  (lib) 

Equation*  (11)  exactly  correspond  to  the  diffusion  aquation  with  an 
instantanaout  point  source  located  at  ($',x')  for  Dlrchlet  and  Neumann 
conditions,  respectively9.  The  solutions  are: 


UT  -  l  2  sln(ntx)sin(ftfx')  e 
*  n-1 


Gj  •  1+  [  2  cos(nix)cos(n*x')  e 

n-i 


B  2  2 
yn-  Pjn* 


ntaCC-C') 


*<*-*'> 


(12a) 


(12b) 


(13) 


The  solution  of  eqs.  (7a)  acd  (7b)  is  now  expressed  in  terms  of  the  Green's 
functions  GT  and  G?  in  the  form  of  integral  equations  containing  the  influence 
of  the  source  terms  and  boundary  and  inlet  conditions: 


-  /J  WVI*—*'  ♦  P1  /o  ''hO’1*'-44' 

"  pi  ll  Jt 

-  /» -  r  n 

where  Che  source  terms  and  Sy  are  given  by: 
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(14b) 
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Solution  of  equations  (14a)  and  (14b)  for  tha  gaa  temperature  and  fuel 
concent  rat  ion,  and  of  equations  (7c)  and  (7d)  for  the  droplet  radius  and 
temperature  requires  an  iterative  scheme  because  the  unknown  functions  appear 
in  the  integrals  of  the  source  terms.  A  more  efficient  numerical  solutioo  can 
be  achieved  by  differentiating  equations  (14a)  and  (14b)  with  respect  to  g,  a 
procedure  that  transforms  then  into  ordinary  differential  equations  which  can 
be  readily  solved  with  a  Runge-Kutta  method.  We  note  the  series-character  of 
the  Green's  functions,  and  therefore  of  the  kernels  in  the  Integral  source 
terns,  and  let  the  functiona  T'  and  Yp  be  composed  of  an  infinite  number  of 
functions  plus  some  known  functiona  related  to  the  inlet  and  boundary 
conditions: 


T»  -  f  ♦  X  2T'nsin(nwx) 
n-1 


(16a) 


’»  *  V  ’to  *  P 


(16b) 


Here.fp  corresponds  to  the  first  three  terms  of  equation  (14a)  while  fp 
corresponds  to  the  first  tern  of  equation  (14b). 

Upon  differentiating  equations  (14a)  and  (14b)  with  respect  to  g  and 
comparing  terms  In  the  Infinite  series  of  equations  (12)  and  (16),  and  after 
some  manipulation  one  obtains: 


ST  +  YnT’n  “  ^I(ST8ln(rhrX,>!x.-Xj 
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(17a) 
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(17b) 


subject  to  th*  conditions: 


T;  -  Tfn  -  0  *  «  “  0 


(17c) 


Equations  (17)  arc  aolvsd  togathar  with  aquatints  (7c)  and  (7d)  using  a  Runga- 
Kutta-Glll  4th  order  scheme  up  to  the  value  of  n  for  which  certain  convergence 
criteria  are  net . 


Results 

The  results  of  the  analysis  will  be  presented  In  the  following  manner. 
First,  we  show  and  discuss  the  results  of  a  characteristic  case,  including 
droplet  temperature  and  else  histories  as  well  as  the  development  of  the  gas- 
phase  temperature  and  fuel-concentrctlon  profiles.  We  continue  the 
presentation  of  the  results  with  an  analysis  of  the  affects  of  the  selected 
parameters,  Pp  P2,  P3,  and  PA. 


The  standard  case  corresponds  to  liquid  n-heptane  droplets  at  300 
K  and  150  pm  in  diameter  injected  with  a  velocity  of  1  m/s  In  an  air  stream 
with  a  linear  temperature  profile  which  goes  from  2050  K  on  the  left  boundary 
to  1000  K  on  the  right  boundary  and  a  velocity  of  5  m/s.  The  dimensionless 
parameters  have  the  following  values:  Pj  »  0.1,  ?2  •  0.1,  P3  •  2,  -  12. 

The  droplet  sheets  are  located  at  one  third  and  two  thirds  of  the  distance 
from  the  left  wall  to  the  right  wall. 

Figures  1  and  2  show  the  droplet-temperature  and  droplet-radius  histories 
on  the  two  droplet  sheets.  The  first  sheet,  1,  is  closer  to  the  hot  boundary 
and  therefore  these  droplets  heat  up  and  vaporize  faster.  There  is  a 
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transient  heating  period  of  roughly  3 3Z  of  Cha  droplet's  llfetlae  during  which 
alaoat  401  of  tha  droplet  la  vaporized  and  after  which  an  alaoat  constant 
droplet  teaperature  la  aatabliehad.  Thia  teaperatura  ia  appreciably  lower 
than  the  aaturatioa  teaperature  at  the  prevailing  praaaura  (T‘  t  ■  0.22). 

Note  that  with  leaa  volatile  fuels,  the  tranalent  behavior  would  be  even  wore 
pronounced. , 

Figure  3  shows  the  development  of  the  fuel  aass  fraction  profiles  in  the 
gas  phase.  The  two  aaxlaa  correspond  to  the  location  of  the  droplet  aheets. 

The  sheet  locked  at  g  *  1/3  is  exposed  to  a  higher  gas  teaperature  and 
therefore  vaporizes  faster  and  produces  aore  fuel  vapor  than  the  sheet  located 
at  g  •  2/3.  However,  after  the  droplets  on  the  first  sheet  are  coapletely 
vaporized  at  5  -  0.725  and  due  to  tha  effect  of  gas  diffusion,  the  fuel 
concentration  at  the  location  of  the  second  sheet  is  evsntually  higher  than 
anywhere  else.  Eventually,  as  g  goes  to  infinity,  a  unifora  profile  results. 

The  gas-phase  teaperature  is  given  in  Figure  4  at  three  different  axial 
locations.  At  the  inlet,  g  •  0,  the  teaperature  profile  is  linear  but  due  to 
heat  transfer  to  the  droplets,  two  energy  sinks  are  developed  at  the  droplet- 
sheet  locations  which  produce  the  two  teaperature  depressions 
near  g  •  1/3  and  g  -  2/3  (g  -  0.2).  Eventually  (g  -  1)  the  two  depressions 
are  smoothed  out  by  the  gas  diffusion  and  as  g  ♦  •,  the  teaperature  profile 
recovers  its  linear  shape. 

Effect  of  Pj 

Figures  5  and  6  correspond  to  a  four-fold  increase  of  P^  (P^  •  0.4)  while 
other  paraaeters  are  kept  unchanged. 

Increasing  Pj  results  in  a  reduction  of  the  tlae  required  by  the  gas- 
diffusion  aechanisa  to  smooth  out  the  differences  created  by  the  heat  losses 
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and  mii  gain*  occurring  at  tha  droplet-sheet  locations.  This  rasults  in  less 
pronounced  aaxiaa  in  the  fuel  aass  fraction  profiles  and  a  smoother 
distribution  of  fuel  vapor  over  the  cross  sectional  area  aa  shorn  in  Figure 
5*  The  droplet-sixe  histories  on  the  two  droplet  sheets  for  this  case  are 
shown  in  Figure  6.  Since  the  droplets  are  now  exposed  to  a  locally  higher 
temperature,  the  vaporisation  tiaes  are  reduced  by  approxlaately  72  with 
respect  to  those  of  Figure  2. 

Effect  of  ?2 

The  Influence  of  Y^  ls  discussed  with  the  aid  of  Figure  7.  This  figure 
corresponds  to  a  four-fold  reduction  of  Pj  (?£  ■  0.025)  while  keeping  all 
other  paraaeters  as  in  the  base  case.  As  explained  in  the  Analysis  section, 

?2  is  the  ratio  of  the  incoming  liquid  and  gas  aass  flows.  Figure  7  shows 
that  a  reduction  in  P2  by  a  factor  of  4  nearly  directly  correlates  with  a 
reduction  by  the  seas  factor  of  the  fuel  vapor  present  in  the  gas  phase. 
Another  consequence  not  presented  here  is  a  reduction  in  the  droplet  lifetiaes 
of  roughly  102  with  respect  to  th*  base  case. 

Effect  of  Pj 

A  reduction  by  a  factor  of  two  of  the  transfer  number  P^  produces  a 
marked  increase  in  the  droplet  lifetime  (452  in  this  case)  as  Indicated  by 
Figure  8  and  a  reduction  of  the  final  droplet  temperature.  This  is  a 
consequence  of  an  increase  in  the  ratio  of  the  heat  required  to  vaporise  the 
droplets  to  the  heat  available  in  the  gas-phase.  The  development  of  the  fuel 
aass  fraction  profile  in  the  gas  phase  takes  longer  here  than  in  the  base  c«*e 
as  can  be  seen  in  Figure  9. 


The  Reynolds  nuaber  appears  In  cha  problaa  through  tha  ralative-valocity 
cor  race  ion  In  tha  aass  vaporization  rata  of  Eqn  (2).  With  tha  aaauaptiona  of 
constant  relative  valocity  between  tha  gaa  and  liquid  phases  and  constant 
physical  properties,  tha  Raynolds  nuabar,  P^,  varies  throughout  the  flow  only 
due  to  variations  in  droplat  size.  Decreasing  the  initial  Reynolds  nuaber 
reduces  tha  aass  vaporization  rata  and  tharafora  lncrsasaa  the  droplet 
llfetlae.  Figures  2  and  10  show  an  Increase  of  droplet  llfetlaes  of 
approximately  25X  when  the  Raynolds  nuabar  is  raducad  froa  12  to  4. 
Accordingly,  coaparlng  Figuras  3  and  11,  it  is  observed  that  the  aaount  of 
fuel  present  in  the  gas  phase  in  tha  early  stages  of  tha  process  (g  *0.2)  is 
about  30%  less  in  tha  lover  Raynolds  nuaber  case.  Note,  however,  that  towards 
the  end  of  tha  droplat  llfetlaes  ({  -  0.8)  that  difference  has  alaost 
disappeared. 

Conclusions 

A  slaplified  but  qualitatively  reasonable  aodel  of  heat  transfer  and 
rapid  vaporization  in  a  two-dimensional  spray  has  been  presented.  Based  on 
tha  assumption  of  unifora  flow  velocities  and  the  neglect  of  viscous  exchange 
of  aoaentua  between  the  gas  and  liquid  phases,  an  analytically  tractable 
system  of  partial  differential  equations  has  been  obtained  which  has  been 
further  transforaed  to  a  system  of  Integral  equations  with  the  use  of  Green's 
functions.  While  the  model  Introduces  considerable  slaplif icatlon  by  removing 
the  nonlinearities  froa  the  differential  operators,  It  retains  the  complexity 
and  nonlinearity  of  the  source  terms,  making  it  possible  to  consider  the 
Inclusion  in  the  aodel  of  the  chemical-reaction  terms.  The  analysis 


Identifies  several  diaensionless  parameters  thee  characterize  the  problea  end 
the  Influence  of  four  of  these  is  presented  In  detail,  namely:  a 
characteristic  vaporisation-diffusion  tine  ratio,  a  liquid-to-gas  ness  flow 
ratio,  a  Spalding  transfer  nuaber,  and  a  Reynolds  number. 

The  aodel  offers  some  quantitatively  crude  but  qualitatively  very 
inforaative  results  of  the  complex  aechanisa  of  heat  and  aass  transfer 
occuring  prior  to  the  ignition  of  two  parallel  fuel  droplet  streaas,  and  in 
particular  it  serves  as  a  base  for  a  simplified  aodel  of  ignition  of  a  droplet 
streaa  by  another  droplet  streaa,  a  topic  that  ia  presently  under 
investigation. 
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reference  transfer  nun bar  defined  in  (9) 
effective  transfer  nuaber 
liquid  specific  beet 
gee  constant-pressure  specific  beat 
convective  correction,  •  1  ♦  0.3  (Re) 
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mss  diffusion  coefficient 

functions  described  in  the  context  of  Eqns  (16a,b) 
Green's  functions  described  by  Eqns  (12e,b) 
channel  width 

effective  beat  of  vaporisation 
latent  heat  of  vaporisation 
droplet  mss 

mss  vaporisation  rata  given  in  (2) 
gas  aolecular  weight 
liquid  Mlecular  weight 
droplet-sheet  density;  series  index 
pressure 

novael  ateospherlc  pressure 
paraaoters  defined  in  (9) 
droplet  radius 

nonulised  droplet  radius,  R'  -  R/R^ 
initial  droplet  radius 
universal  gas  constant 

2RI  V°t  I 

Reynolds  nuaber,  Ra  •  - - — *■ 

v 

source  teras  described  by  Eqns  (15a, b) 
gas  teaperature 
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normalized  gas  temperature,  T1  •  (T-T  ^/(T^T^) 

inlet  gaa  temperature 

cold  wall  taeparaturc 

hot  wall  tsaparacura 

liquid  taaparatura 

normalized  liquid  taaparatura,  T*  ■  (I  -I  . )/(T  -T  ) 

t  I  i*  •  •• 

n-th  coaponent  of  taaparatura  sariaa  (Eq  16a) 
liquid  velocity 
tranavarse  coordinate 
fuel  aaaa  fraction 

fuel  aaaa  fraction  at  droplat  surface 

n-th  coaponent  of  fuel  aaaa  fraction  sariaa  (Eq  16b) 

atreaawlsa  coordinate 


Greek  ayabols 
a 

*n 

6 


thernal  dlffuaivlty 
paraaatar  defined  in  (13) 

Dirac  delta  function  $  ■  j(x-x^) 

diaenaionless  Dirac  delta  function  6’  -  j'(x“Xj)  *  H 
klneaatic  viscosity 

diaenaionless  streaawlse  coordinate  (Eqn  6a) 
gas  density 
liquid  density 

characteristic  vaporization  time  defined  in  (6a) 
noraalized  transverse  coordinate,  y  ■  x/L 
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Subscripts 
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Figure  Cape ions 

1.  Normalized  droplet  temperature  histories  for  the  hese  ease.  I:  Droplet 
sheet  near  hot  vail,  11:  Droplet  sheet  near  cold  wall 

2.  Normalized  droplet  else  histories  for  the  base  case.  I:  Droplet  sheet 
near  hot  wall,  II:  Droplet  sheet  near  cold  wall 

3.  Development  of  the  gas-phase  fuel  mass  fraction  profiles  for  the  base  case 

4.  Development  of  the  gas-phase  temperature  profiles  for  the  base  case 

3.  Development  of  the  gas-phase  fuel  mass  fraction  profiles  for  the  case  of 
*1  -  0.4 

6.  Normalized  droplet  size  histories  for  the  case  of  Pj  »  0.4 

7.  Development  of  the  gas-phase  fuel  mass  fraction  profiles  for  the  case  of 
P2  -  0.025 

8.  Normalized  droplet  size  histories  for  the  case  of  Pj  -  1 

9.  Development  of  the  gas-phase  fuel  mass  fraction  profllee  for  the  case  of 

P3  -  1 

10.  Noraalized  droplet  size  histories  for  the  case  of  P^  -  4 

11.  Development  of  the  gas-phase  fuel  mass  fraction  profiles  for  the  case  of 
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